Analysis of ozone and nitric acid in spring and summer Arctic pollution using aircraft, ground-based, satellite observations and MOZART-4 model: source attribution and partitioning

In this paper, we analyze tropospheric O_3 together with HNO_3 during the POLARCAT (Polar Study using Aircraft, Remote Sensing, Surface Measurements and Models, of Climate, Chemistry, Aerosols, and Transport) program, combining observations and model results. Aircraft observations from the NASA ARCTAS (Arctic Research of the Composition of the Troposphere from Aircraft and Satellites) and NOAA ARCPAC (Aerosol, Radiation and Cloud Processes affecting Arctic Climate) campaigns during spring and summer of 2008 are used together with the Model for Ozone and Related Chemical Tracers, version 4 (MOZART-4) to assist in the interpretation of the observations in terms of the source attribution and transport of O_3 and HNO_3 into the Arctic (north of 60° N). The MOZART-4 simulations reproduce the aircraft observations generally well (within 15%), but some discrepancies in the model are identified and discussed. The observed correlation of O_3 with HNO_3 is exploited to evaluate the MOZART-4 model performance for different air mass types (fresh plumes, free troposphere and stratospheric-contaminated air masses). Based on model simulations of O_3 and HNO_3 tagged by source type and region, we find that the anthropogenic pollution from the Northern Hemisphere is the dominant source of O3 and HNO3 in the Arctic at pressures greater than 400 hPa, and that the stratospheric influence is the principal contribution at pressures less 400 hPa. During the summer, intense Russian fire emissions contribute some amount to the tropospheric columns of both gases over the American sector of the Arctic. North American fire emissions (California and Canada) also show an important impact on tropospheric ozone in the Arctic boundary layer. Additional analysis of tropospheric O_3 measurements from ground-based FTIR and from the IASI satellite sounder made at the Eureka (Canada) and Thule (Greenland) polar sites during POLARCAT has been performed using the tagged contributions. It demonstrates the capability of these instruments for observing pollution at northern high latitudes. Differences between contributions from the sources to the tropospheric columns as measured by FTIR and IASI are discussed in terms of vertical sensitivity associated with these instruments. The first analysis of O_3 tropospheric columns observed by the IASI satellite instrument over the Arctic is also provided. Despite its limited vertical sensitivity in the lowermost atmospheric layers, we demonstrate that IASI is capable of detecting low-altitude pollution transported into the Arctic with some limitations

Atmospheric transport and chemistry of trace gases in LMDz5B: evaluation and implications for inverse modelling

Representation of atmospheric transport is a major source of error in the estimation of greenhouse gas sources and sinks by inverse modelling. Here we assess the impact on trace gas mole fractions of the new physical parameterizations recently implemented in the atmospheric global climate model LMDz to improve vertical diffusion, mesoscale mixing by thermal plumes in the planetary boundary layer (PBL), and deep convection in the troposphere. At the same time, the horizontal and vertical resolution of the model used in the inverse system has been increased. The aim of this paper is to evaluate the impact of these developments on the representation of trace gas transport and chemistry, and to anticipate the implications for inversions of greenhouse gas emissions using such an updated model. Comparison of a one-dimensional version of LMDz with large eddy simulations shows that the thermal scheme simulates shallow convective tracer transport in the PBL over land very efficiently, and much better than previous versions of the model. This result is confirmed in three-dimensional simulations, by a much improved reproduction of the radon-222 diurnal cycle. However, the enhanced dynamics of tracer concentrations induces a stronger sensitivity of the new LMDz configuration to external meteorological forcings. At larger scales, the inter-hemispheric exchange is slightly slower when using the new version of the model, bringing them closer to observations. The increase in the vertical resolution (from 19 to 39 layers) significantly improves the representation of stratosphere/troposphere exchange. Furthermore, changes in atmospheric thermodynamic variables, such as temperature, due to changes in the PBL mixing modify chemical reaction rates, which perturb chemical equilibriums of reactive trace gases. One implication of LMDz model developments for future inversions of greenhouse gas emissions is the ability of the updated system to assimilate a larger amount of high-frequency data sampled at high-variability stations. Others implications are discussed at the end of the paper

Technical Note: Ozonesonde climatology between 1995 and 2011: description, evaluation and applications

An ozone climatology based on ozonesonde measurements taken over the last 17 yr has been constructed for model evaluation and comparisons to other observations. Vertical ozone profiles for 42 stations around the globe have been compiled for the period 1995–2011, in pressure and tropopause-referenced altitudes. For each profile, the mean, standard deviation, median, the half-width are provided, as well as information about interannual variability. Regional aggregates are formed in combining stations with similar ozone characteristics. The Hellinger distance is introduced as a new diagnostic to identify stations that describe similar shapes of ozone probability distribution functions (PDFs). In this way, 12 regions were selected covering at least 2 stations and the variability among those stations is discussed. Significant variability with longitude of ozone distributions in the troposphere and lower stratosphere in the northern mid- and high latitudes is found. The representativeness of regional aggregates is discussed for high northern latitudes, Western Europe, Eastern US, and Japan, using independent observations from surface stations and MOZAIC aircraft data. Good agreement exists between ozonesondes and aircraft observations in the mid-troposphere and between ozonesondes and surface observations for Western Europe. For Eastern US and high northern latitudes, surface ozone values from ozonesondes are biased 10 ppb high compared to independent measurements. An application of the climatology is presented using the NCAR CAM-Chem model. The climatology allows evaluation of the model performance regarding ozone averages, seasonality, interannual variability, and the shape of ozone distributions. The new assessment of the key features of ozone distributions gives deeper insights into the performance of models

Inter-model comparison of global hydroxyl radical (OH) distributions and their impact on atmospheric methane over the 2000–2016 period

The modeling study presented here aims to estimate how uncertainties in global hydroxyl radical (OH) distributions, variability, and trends may contribute to resolving discrepancies between simulated and observed methane (CH4) changes since 2000. A multi-model ensemble of 14 OH fields was analyzed and aggregated into 64 scenarios to force the offline atmospheric chemistry transport model LMDz (Laboratoire de Meteorologie Dynamique) with a standard CH4 emission scenario over the period 2000–2016. The multi-model simulated global volume-weighted tropospheric mean OH concentration ([OH]) averaged over 2000–2010 ranges between 8:7*10^5 and 12:8*10^5 molec cm-3. The inter-model differences in tropospheric OH burden and vertical distributions are mainly determined by the differences in the nitrogen oxide (NO) distributions, while the spatial discrepancies between OH fields are mostly due to differences in natural emissions and volatile organic compound (VOC) chemistry. From 2000 to 2010, most simulated OH fields show an increase of 0.1–0:3*10^5 molec cm-3 in the tropospheric mean [OH], with year-to-year variations much smaller than during the historical period 1960–2000. Once ingested into the LMDz model, these OH changes translated into a 5 to 15 ppbv reduction in the CH4 mixing ratio in 2010, which represents 7%–20% of the model-simulated CH4 increase due to surface emissions. Between 2010 and 2016, the ensemble of simulations showed that OH changes could lead to a CH4 mixing ratio uncertainty of > 30 ppbv. Over the full 2000–2016 time period, using a common stateof- the-art but nonoptimized emission scenario, the impact of [OH] changes tested here can explain up to 54% of the gap between model simulations and observations. This result emphasizes the importance of better representing OH abundance and variations in CH4 forward simulations and emission optimizations performed by atmospheric inversions

Impact of transport model errors on the global and regional methane emissions estimated by inverse modelling

A modelling experiment has been conceived to assess the impact of transport model errors on methane emissions estimated in an atmospheric inversion system. Synthetic methane observations, obtained from 10 different model outputs from the international TransCom-CH[subscript 4] model inter-comparison exercise, are combined with a prior scenario of methane emissions and sinks, and integrated into the three-component PYVAR-LMDZ-SACS (PYthon VARiational-Laboratoire de Météorologie Dynamique model with Zooming capability-Simplified Atmospheric Chemistry System) inversion system to produce 10 different methane emission estimates at the global scale for the year 2005. The same methane sinks, emissions and initial conditions have been applied to produce the 10 synthetic observation datasets. The same inversion set-up (statistical errors, prior emissions, inverse procedure) is then applied to derive flux estimates by inverse modelling. Consequently, only differences in the modelling of atmospheric transport may cause differences in the estimated fluxes. In our framework, we show that transport model errors lead to a discrepancy of 27 Tg yr[superscript −1] at the global scale, representing 5% of total methane emissions. At continental and annual scales, transport model errors are proportionally larger than at the global scale, with errors ranging from 36 Tg yr[superscript −1] in North America to 7 Tg yr[superscript −1] in Boreal Eurasia (from 23 to 48%, respectively). At the model grid-scale, the spread of inverse estimates can reach 150% of the prior flux. Therefore, transport model errors contribute significantly to overall uncertainties in emission estimates by inverse modelling, especially when small spatial scales are examined. Sensitivity tests have been carried out to estimate the impact of the measurement network and the advantage of higher horizontal resolution in transport models. The large differences found between methane flux estimates inferred in these different configurations highly question the consistency of transport model errors in current inverse systems. Future inversions should include more accurately prescribed observation covariances matrices in order to limit the impact of transport model errors on estimated methane fluxes

Reconciling the bottom-up and top-down estimates of the methane chemical sink using multiple observations

The methane chemical sink estimated by atmospheric chemistry models (bottom-up method) is significantly larger than estimates based on methyl chloroform (MCF) inversions (top-down method). The difference is partly attributable to large uncertainties in hydroxyl radical (OH) concentrations simulated by the atmospheric chemistry models used to derive the bottom-up estimates. In this study, we propose a new approach based on OH precursor observations and a chemical box model. This approach contributes to improving the 3D distributions of tropospheric OH radicals obtained from atmospheric chemistry models and reconciling bottom-up and top-down estimates of the chemical loss of atmospheric methane. By constraining simulated OH precursors with observations, the global mean tropospheric column-averaged air-mass-weighted OH concentration ([OH]trop-M) is ∼10×105 molec. cm−3 (which is 2×105 molec. cm−3 lower than the original model-simulated global [OH]trop-M) and agrees with that obtained by the top-down method based on MCF inversions. With OH constrained by precursor observations, the methane chemical loss is 471–508 Tg yr−1, averaged from 2000 to 2009. The new adjusted estimate is in the range of the latest top-down estimate of the Global Carbon Project (GCP) (459–516 Tg yr−1), contrary to the bottom-up estimates that use the original model-simulated OH fields (577–612 Tg yr−1). The overestimation of global [OH]trop-M and methane chemical loss simulated by the atmospheric chemistry models is caused primarily by the models' underestimation of carbon monoxide and total ozone column, and overestimation of nitrogen dioxide. Our results highlight that constraining the model-simulated OH fields with available OH precursor observations can help improve bottom-up estimates of the global methane sink.</p

On the use of Earth Observation to support estimates of national greenhouse gas emissions and sinks for the Global stocktake process: lessons learned from ESA-CCI RECCAP2

The Global Stocktake (GST), implemented by the Paris Agreement, requires rapid developments in the capabilities to quantify annual greenhouse gas (GHG) emissions and removals consistently from the global to the national scale and improvements to national GHG inventories. In particular, new capabilities are needed for accurate attribution of sources and sinks and their trends to natural and anthropogenic processes. On the one hand, this is still a major challenge as national GHG inventories follow globally harmonized methodologies based on the guidelines established by the Intergovernmental Panel on Climate Change, but these can be implemented differently for individual countries. Moreover, in many countries the capability to systematically produce detailed and annually updated GHG inventories is still lacking. On the other hand, spatially-explicit datasets quantifying sources and sinks of carbon dioxide, methane and nitrous oxide emissions from Earth Observations (EO) are still limited by many sources of uncertainty. While national GHG inventories follow diverse methodologies depending on the availability of activity data in the different countries, the proposed comparison with EO-based estimates can help improve our understanding of the comparability of the estimates published by the different countries. Indeed, EO networks and satellite platforms have seen a massive expansion in the past decade, now covering a wide range of essential climate variables and offering high potential to improve the quantification of global and regional GHG budgets and advance process understanding. Yet, there is no EO data that quantifies greenhouse gas fluxes directly, rather there are observations of variables or proxies that can be transformed into fluxes using models. Here, we report results and lessons from the ESA-CCI RECCAP2 project, whose goal was to engage with National Inventory Agencies to improve understanding about the methods used by each community to estimate sources and sinks of GHGs and to evaluate the potential for satellite and in-situ EO to improve national GHG estimates. Based on this dialogue and recent studies, we discuss the potential of EO approaches to provide estimates of GHG budgets that can be compared with those of national GHG inventories. We outline a roadmap for implementation of an EO carbon-monitoring program that can contribute to the Paris Agreement

Evaluating methane inventories by isotopic analysis in the London region

A thorough understanding of methane sources is necessary to accomplish methane reduction targets. Urban environments, where a large variety of methane sources coexist, are one of the most complex areas to investigate. Methane sources are characterised by specific δ13C-CH4 signatures, so high precision stable isotope analysis of atmospheric methane can be used to give a better understanding of urban sources and their partition in a source mix. Diurnal measurements of methane and carbon dioxide mole fraction, and isotopic values at King’s College London, enabled assessment of the isotopic signal of the source mix in central London. Surveys with a mobile measurement system in the London region were also carried out for detection of methane plumes at near ground level, in order to evaluate the spatial allocation of sources suggested by the inventories. The measured isotopic signal in central London (−45.7 ±0.5‰) was more than 2‰ higher than the isotopic value calculated using emission inventories and updated δ13C-CH4 signatures. Besides, during the mobile surveys, many gas leaks were identified that are not included in the inventories. This suggests that a revision of the source distribution given by the emission inventories is needed

Estimating methane emissions in the Arctic nations using surface observations from 2008 to 2019

The Arctic is a critical region in terms of global warming. Environmental changes are already progressing steadily in high northern latitudes, whereby, among other effects, a high potential for enhanced methane (CH4) emissions is induced. With CH4 being a potent greenhouse gas, additional emissions from Arctic regions may intensify global warming in the future through positive feedback. Various natural and anthropogenic sources are currently contributing to the Arctic's CH4 budget; however, the quantification of those emissions remains challenging. Assessing the amount of CH4 emissions in the Arctic and their contribution to the global budget still remains challenging. On the one hand, this is due to the difficulties in carrying out accurate measurements in such remote areas. Besides, large variations in the spatial distribution of methane sources and a poor understanding of the effects of ongoing changes in carbon decomposition, vegetation and hydrology also complicate the assessment. Therefore, the aim of this work is to reduce uncertainties in current bottom-up estimates of CH4 emissions as well as soil oxidation by implementing an inverse modelling approach in order to better quantify CH4 sources and sinks for the most recent years (2008 to 2019). More precisely, the objective is to detect occurring trends in the CH4 emissions and potential changes in seasonal emission patterns. The implementation of the inversion included footprint simulations obtained with the atmospheric transport model FLEXPART (FLEXible PARTicle dispersion model), various emission estimates from inventories and land surface models, and data on atmospheric CH4 concentrations from 41 surface observation sites in the Arctic nations. The results of the inversion showed that the majority of the CH4 sources currently present in high northern latitudes are poorly constrained by the existing observation network. Therefore, conclusions on trends and changes in the seasonal cycle could not be obtained for the corresponding CH4 sectors. Only CH4 fluxes from wetlands are adequately constrained, predominantly in North America. Within the period under study, wetland emissions show a slight negative trend in North America and a slight positive trend in East Eurasia. Overall, the estimated CH4 emissions are lower compared to the bottom-up estimates but higher than similar results from global inversions.</p

Inverse modelling of European CH4 emissions during 2006-2012 using different inverse models and reassessed atmospheric observations

We present inverse modelling (top down) estimates of European methane (CH4) emissions for 2006-2012 based on a new quality-controlled and harmonised in situ data set from 18 European atmospheric monitoring stations. We applied an ensemble of seven inverse models and performed four inversion experiments, investigating the impact of different sets of stations and the use of a priori information on emissions.The inverse models infer total CH4 emissions of 26.8 (20.2-29.7) TgCH(4) yr(-1) (mean, 10th and 90th percentiles from all inversions) for the EU-28 for 2006-2012 from the four inversion experiments. For comparison, total anthropogenic CH4 emissions reported to UNFCCC (bottom up, based on statistical data and emissions factors) amount to only 21.3 TgCH(4) yr(-1) (2006) to 18.8 TgCH(4) yr(-1) (2012). A potential explanation for the higher range of top-down estimates compared to bottom-up inventories could be the contribution from natural sources, such as peatlands, wetlands, and wet soils. Based on seven different wetland inventories from the Wetland and Wetland CH4 Inter-comparison of Models Project (WETCHIMP), total wetland emissions of 4.3 (2.3-8.2) TgCH(4) yr(-1) from the EU-28 are estimated. The hypothesis of significant natural emissions is supported by the finding that several inverse models yield significant seasonal cycles of derived CH4 emissions with maxima in summer, while anthropogenic CH4 emissions are assumed to have much lower seasonal variability. Taking into account the wetland emissions from the WETCHIMP ensemble, the top-down estimates are broadly consistent with the sum of anthropogenic and natural bottom-up inventories. However, the contribution of natural sources and their regional distribution remain rather uncertain.Furthermore, we investigate potential biases in the inverse models by comparison with regular aircraft profiles at four European sites and with vertical profiles obtained during the Infrastructure for Measurement of the European Carbon Cycle (IMECC) aircraft campaign. We present a novel approach to estimate the biases in the derived emissions, based on the comparison of simulated and measured enhancements of CH4 compared to the background, integrated over the entire boundary layer and over the lower troposphere. The estimated average regional biases range between -40 and 20% at the aircraft profile sites in France, Hungary and Poland.</p